DYNAMICAL SYSTEMS OF SIMPLICES IN DIMENSION 2 OR 3 
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Abstract. Let To = (Aq---Aq) be a d-simplex, Go its centroid, S its circumsphere, O the 
center of S. Let (A\) be the points where S intersects the lines (GoA ), T\ the d-simplex 
(Ai---Af), and Gi its centroid. By iterating this construction, a dynamical system of d- 
f-') ■ simplices (T) with centroids (Gi) is constructed. For d = 2 or 3, we prove that the sequence 

' (OGi)i is decreasing and tends to 0. We consider the sequences (Tn)i and {Ta+i)i\ for d = 2 

£Nj they converge to two equilateral triangles with at least quadratic speed; for d — 3 they converge 

to two isosceles tetrahedra with at least geometric speed. In this last case, we give an explicit 
expression of the lengths of the edges of the limit form. We show also that if To is a planar cyclic 
quadrilateral then (T n ) converges to a rectangle with at least geometric speed or eventually to 
a square with a speed that is conjectured as cubic. The proofs are largely algebraic and use 
Grobner basis computations. 
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1. Introduction 



Q 

^^^^ 

1.1. The general problem. Let Tq = (Aq---Aq) be a d-simplex, Go its centroid, S its cir- 
cumsphere in M d , O the center of S. Let (A\) be the points where S intersects the lines (GoA ) 
respectively, T\ the d-simplex (A® ■ ■ ■ Af), and G\ its centroid. If we iterate this construction 
then we produce a dynamical system of ci-simplices (%) with centroids (Gi). 

Let M, N be two points of M d ; MN refers to the euclidean norm of the vector MN. 
■ Numerical investigations, using Maple, with thousands of random simplices in dimensions up 

£T) . to 20 indicate that: 

^ , i) the sequence (OGi)i is decreasing and tends to (from ||). 

ii) the sequences (72i)« and (^2i+i)j converge to two d-simplices with centroid O. 

o ; 

Q\ [ Remark. The condition "circumcenter=centroid" for a d-simplex is equivalent to the statement 

that the sum of the squares of the edge lengths of each facet is the same for all d+ 1 facets. 

We prove the assertions i) and ii) in the cases d = 2 and d = 3. We conjecture that the result 
^ . is valid for any d. 

H ' From now on, we assume that the radius of S is 1. 

& '_ 

1.2. The case d = 3. % = (AiBiCiDi) is a tetrahedron such that its four vertices are not 
coplanar. These vertices are ordered and % is viewed as a point of S 4 
Let (f> and tp be the functions which transform Tq into T\ and G\. 

Remark, i) (f> does not admit any non planar tetrahedron as fixed point, 
ii) If T is an isosceles tetrahedron then T is a fixed point of (f>o <j>. 

Our main result is the following: 

The sequences (72i)ieN and (72i+i)iGN are well defined and converge, with at least geometric 
speed, to two non planar isosceles tetrahedra that are symmetric with respect to O. 

We consider also the degenerate case where 7o is a planar cyclic convex quadrilateral; then 
the sequence (%)i converges, with at least geometric speed, to a rectangle. If Tq is a harmonic 
quadrilateral then the limit form is a square and we conjecture that the convergence is with 
order three. 

Moreover, in both cases, we give explicit expressions for the lengths of the edges of the limit 
form from the ones of Tq. 
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1.3. The case d = 2. %,= (AiB{Ci) is a triangle such that its vertices are pairwise distinct. 
We prove that the sequences (T^XeM and (72i+i)iGN are well defined and converge to two 

equilateral triangles that are symmetric with respect to O. Moreover, these sequences converge 
with at least quadratic speed. 

Thus if the dimension is two, it can be observed a much more quick convergence that in the 
dimension three case. 

1.4. Method used. Along this paper, we use rational functions of the square of the lengths 
of the edges of the d-simplices and naturally some systems of polynomial equations appear. 
Computations with such systems requires Grobner basis softwares. For this, we have used the 
J. C. Faugere's software "FGb" (see and the computer algebra system Maple (sec ||) 
which provides some other tools we needed. 



2. Standard definitions and results about tetrahedra. 



2.1. General tetrahedra. Let T = (ABCD) be a tetrahedron, G its centroid and V(T) its 
volume. We assume A,B,C,D are not coplanar and the circumsphere of T has center O and 
radius 1. Let BC = a,CA = b, AB = c, AD = a', BD = b', CD = d. 

Remark that the 3-tuple (a,b,c) does not play the same role as (a',b',c r ). 

Proposition 1. For all T, we have 
1 



i) OG 1 



16 
1 



[a? + b 2 + c 2 + a' 2 + b' 2 + c'' 



(GA 2 + GB 2 + GC 2 + GD 2 ). (See @ p. 64 and ii)). 



ii) GA 2 = 


3 


16 


GB 2 = 


3 




16 


GC 2 = 


3 




16 


GD 2 = 


3 




16 


If U, V, w 





+ c 2 + a 2 ) 



5« 
1 

1 



J2\ 



9 + c' 2 + a' 2 ) 



c' 2 +a 2 + b 2 )-—(c 2 +a ,2 + b' 2 ), 
16 

y 2 + 6 /2 +c /2 ) __L (a 2 + 6 2 + c 2 ) _ 



see 



11 1 p. xv). 



then Gram(u, v, w) refers to the Gram matrix of these three vectors. If U is 
a square matrix then det(U) refers to its determinant. 



Proposition 2. Let T(A, B,C, D) = det 1 c 2 a 2 b' 2 be the Cayley-Menger deter- 
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minant, and let A(A, B, C, D) = det 



With these notations, we have 

i) The volume as a function of a,b,c,a' ,b' ,c' (see ||] p. 168): 

288 x V(T) 2 = 8 x det(Gram(AB, AC, AD)) = F(A, B, C, D). 
Therefore, A,B,C,D are not coplanar if and only ifT{A,B,C,D) > 0. 

ii) Let R be the circumradius of the tetrahedron (here R= 1). 

a) Crelle's formula (1821): 6RV(T) = 5 where 5 denotes the area of the triangle with sides 
of lengths aa',bb',cc' (see [g] p. 250) 
?2 _ A(A,B,C,D) 



b) R 2 



2T(A,B,C,D) 



(see M P- 168). 
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iii) There exists 


a relation between a, b, c, a',b', d : det 
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c 2 a 2 


b' 2 
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b 2 a 2 


c' 2 






^ 1 


a' 2 b' 2 c' 2 


°/ 


Proof, iii) comes 


from ii)b) and R = 1: T(A, B, C, D) + -A(A, B, C, D) = 


0. 



0. 



□ 

Proposition 3. We assume T is a non planar tetrahedron or a cyclic quadrilateral. 
The expression Pt(T) = —A(A,B,C,D) has the following properties: 

i) Pt(T) = (bb' + cc' - aa'){cc' + aa' - bb'){aa' + bb' - cc')(aa' + bb' + cc') 

= 2a 2 a ,2 b 2 b' 2 + 2b 2 b' 2 c 2 c' 2 + 2c 2 c' 2 a 2 a' 2 - aV 4 - 6 4 6' 4 - c 4 c' 4 . 

ii) If T is not planar then each factor of the product appearing in i) is positive (Ptolemy's 
inequality, see (see also pj p. 549, 555). 

Moreover Pt{T) = 576 x V{Ty x. R 2 (From Proposition || 

iii) If T is a cyclic quadrilateral then Pt{T) = 0. 

Conversely let a,b,c,a',b',c' be six positive reals. Does there exist a tetrahedron such that 
the lengths of its edges are these reals? 

Proposition 4. The necessary and sufficient condition that the six lengths form a tetrahedron 
( eventually planar) appears to be: 



i) a,b,c are the lengths of a triangle that is 16 x S 2 
(then S is the area of this triangle). 
/0 1 1 1 1 \ 



2a 2 b 2 + 2b 2 c 2 + 2c 2 a 2 - a 4 - 6 4 - c 4 > 



ii) det 



1 c 2 b 2 a' 2 



1 
1 

V 1 



cf 
b 2 

-J2 







a- b' 2 c' z J 

Proof. (See also [10] p. 545) Assertion i) of Proposition || implies that the condition ii) is 
necessary. 

The following can be generalized in any dimension: 

The tetrahedron exists if and only if Gram(AB, AC, AD) is a symmetric non negative matrix. 
That gives the conditions i),ii). □ 

2.2. Isosceles (or equifacetal) tetrahedra. 

Definition. Let T be a non planar tetrahedron whose circumradius is 1. The tetrahedron T is 
said to be isosceles (or equifacetal) if a = a',b = b',c = c' . 

Proposition 5. We have 

i) T is isosceles if and only if G = O. (see [ffl] p. 197). 

ii) If T is isosceles then 72 V 2 (T) = (b 2 + c 2 - a 2 )(c 2 + a 2 
b 2 < c 2 + a 2 , c 2 <a 2 + b 2 . (see [§ p. 205). 

iii) For all T , < aa' + bb' + cc' < 8. Moreover aa' + bb' + cc' 
(see Q p. 558). 

The following result is a direct consequence of propositions H and |]. 

Proposition 6. Let T = (ABCD) be a tetrahedron eventually planar whose centroid is G. 
GA = GB = GC = GD if and only if Tis an isosceles tetrahedron or a planar rectangle. 

3. Deformation from T to 71 (d = 3) 

3.1. Parameters and notations. We adapt our preceding notations; Aq, Bq, Cq, Dq are four 
points of S. We assume only Aq, Bq, Cq, Dq are not all equal. Then the functions <j>,ij) are 
continuous in Tq. 
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a 


J2 



b' 2 

J2 



> 



b 2 ){a 2 + b 2 - c 2 ) and a 2 < b 2 + c 2 , 
8 if and only ifT is isosceles. 



Remark, i) The following calculations are valid if Aq, Bq, Cq, Dq are four points, not all 

equal, of C, a planar circle of center O and radius 1. 

ii) Let A = {(P, P,P,P) : P £ S}. <j> is defined on 5 4 \A. 

Definition. The parameters (dij)i < j of Tq are the square of the lengths of their edges: d% 2 = 
AoBo 2 ,di3 = AoCo 2 ,di4 = AoDo 2 ,d 23 = B(jCo 2 ,d 2 4 = BqDq 2 , d 3 4 = CqDq 2 . Let us recall that 
they are linked by a polynomial equality (see Assertion Hi) of Proposition [|j and two polynomial 
inequalities (see Proposition^). 

Let gi = G A 2 ,g 2 = G B 2 ,g 3 = G C 2 ,g4 = G D 2 and let p = 1 - OG 2 = -{g x + g 2 + 

93 + ff4) an d Pi be the opposites of the powers of Go and G\ with respect to S. According to 
Proposition 1, (gi)% and po are polynomial functions of (dij)ij. 



Proposition 7. We have 

i) The parameters of T\ are id], = Po 2 — — | 

V 3 <l,9j) :i 



■■) Po 2 ^ dij 



16 ^ l<3 g igj 

Hi) OGq 2 and OG\ 2 are rational functions of the {dij)ij. 



iv) V(Tq) 2 andViTi) 2 are rational functions of the (di 3ll3 . 
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Proof For i): G A x G A 1 = p and G ^i 2 = — . The triangles (G A B ) and (GqBiAi) 

91 

are similar; then — - — - = — - — - and A\B 2 = d\ 2 ^ 



AqB GqBq gig 2 
We deduce ii) from i) and Hi) from ii). We deduce iv) from i) and Proposition g. i). □ 

3.2. Inequalities. 

Let T be a tetrahedron or a cyclic quadrilateral. 

Let "P(T) be the property: "the edges of T satisfy d\ 2 = d 3 ±, d\ 3 = d 2 4, du = d 23 " i.e. "T is 
an isosceles tetrahedron or a planar rectangle". 
The first key is the following result: 



Theorem 1. We have the following inequalities: 
i) OG x < OG ; 



ii) for any (ijkl) G {(1234), (1324), (1423)}, d i3 d M < dU\ 



Hi) Pt(%) < Pt(Ti) (see Proposition 3); 
Moreover, fori), ii) and Hi), equalities stands if and only ifViT^o) holds. 

Proof. • For i): the following choice of unknowns enables us to conclude; we may assume 
9i - 92 = 51,52 ~ 93 = s 2 ,g 3 - #4 = s 3 ,giS A = 1 where s 1 ,s 2 ,s 3 > and s 4 > 0. Let's set 
E = 61 (Pi -Po)9i9293 
Po 

We note that signum(E) = signumiOG^ 2 — OG\ 2 ) and we obtain, using the J. C. Faugere's 
software "FGb", this miraculous result: 

E = d 23 (si - S3) 2 + 16s 2 2 S3 2 S4 + 4d 34 sis 2 + 12sis 2 2 + 4sis 2 s 3 + 20s 2 s 3 2 + 21s 2 s 3 3 s 4 + d 2i s 3 2 + 

d 3 4Si 2 +4siS3 2 +d 2 4Si 2 S3S4+d34Si 2 S3S4 + 12s3 3 +2d 2 4SiS 2 S 3 S4+8s 2 3 +3(i 2 4SiS3 2 S4+6d34SiS 
4c?34S 2 2 + d 3 4,S 3 2 + 5d 3 4S 2 2 S 3 S4 + d 2 4,S\ 2 + 3d 34 SiS 2 2 S4 + 4s 2 3 S 3 S 4 + 3(i 3 4S 2 S3 2 S4 + 2d 24 SiS 3 + 

12s 2 2 s 3 + 2d 34 sis 3 + 4d 34 s 2 s 3 + 9s 3 4 s 4 + 2d 34 s 2 3 s 4 + <i34Si 2 s 2 s 4 + si 2 s 2 s 3 s 4 + 4sis 2 2 s 3 s 4 + 
si 2 S3 2 S4 + 10sis 2 S3 2 S4 + 6siS3 3 S4 + 3d34SiS3 2 S4; obviously E is a non negative real. 

Moreover if E = then s\ = s 2 = s 3 = and, according to Proposition ||, V(Tq) holds. 

4 

• For ii): let A = — ^ . By Proposition PtL d)Al, = Adijd kl . According to the AM-GM 

9192939a ^ 3 

inequality, A > 1 and dijdki < djjd^. 



Moreover if d^d^i = djjd^ then A = 1 and according to the properties of the AM-GM 
inequality, gi = g 2 = g 3 = g A and V{%) holds. 

• For Hi): An easy computation gives Pt(T~\) = h?Pt(To) and we reason as for ii). □ 

4. Solution of the case d = 3. Part 1 
The (d^)j<j and ((%)i<j refer to the parameters of T n and Tq. 

Theorem 2. i) The sequence (OG n ) n tends to 0. 
ii) Let T = (ABCD) be a cluster point of the bounded sequence (%). Then T is not flat and 
is isometric to a fixed isosceles tetrahedron that admits the following parameters: 
dyi = Ld\i ^34,^13 2 = Td\zd2A,d\A = ^^14^23 where 
L- 64 

(Vdl2^34 + 7^13^24 + V^14^23) 2 

Proof. According to Theorem [j] Hi), the sequence (Pt{%)) is increasing, then is positive and, 
for all i, Ai, B>i,Ci, Di are not coplanar; for all i, Gi ^ {Ai, Bi, d, Di} then the sequence {%) 
is well defined and the functions <j),ip are continuous in (%). Moreover the bounded sequence 
(Pt(Ti)) converges to Pt > 0. 

According to Theorem [l] i), the bounded sequence (OG n ) is decreasing and converges to r > 0. 
We can extract a subsequence (T nk )k such that (A nk )k, (B nk )y., [C nk )y., (D nk )k converge to 
A, B, C, D. Let T be the tetrahedron (ABCD) and G its centroid. Therefore (G nk )k converges 
to G,OG = r and Pt(T) = Pt; thus A, B, C, D are not coplanar and 4>(T) = T ' = (A'B'C'D 1 ) 
is well defined; let G' be its centroid. 

Assume T is not isosceles; then Theorem [j] i) implies that OG' < OG. Let e E (0, OG — OG') 
and let a > such that if \\T — T n \\ < a then OG n +i — OG' < e. There exists an integer n/% 
such that \\T — T nk \\ < ot; then OG nk+ \ — OG' < e. Thus OG rik+ \ < OG = r, a contradiction. 
Therefore T is isosceles and (OG n ) n converges to 0. A parameter df? of T is the limit of the 
sequence (d" fe ) fc . 

According to the proof of Theorem [l] ii), if {i,j,k,l} is a permutation of {1,2,3,4}, then 
the bounded sequence (d^-^)n satisfy d^^d 1 ^ 1 = A n d^d^; it is increasing then is convergent. 
Thus the infinite product Il^ A n converges to L > 1 such that d^ 2 = Ldijdki- 
^12 + ^13 + ^14 = 8 gives the explicit value of L as a function in the (dij). 
Thus there exists an unique value of d^ that is valid for all cluster points. □ 

Corollary 1. If i < j then the sequence (dfj) n converges to d^ . 

Proof Let U = {(1234), (1324), (1423)}. If (ijkl) £ U then converges to d°° 2 and 

E(iifci)6« \J d ij d kl converges to 8. 

2 E(^)e« ^&%<Ki < EiKjdfj < 16; thus E( ijk l)eU (\[<%j ~ \/dJi) converges to and if 
(ijkl) G then d^- — dfy converges to 0; therefore dfj and converge to d^. □ 

It remains to prove that the (T~2i) cannot turn around O. 

5. Solution of the case d = 3. Part 2 

We assume Tq is not isosceles. We study the convergence speed of the sequence (OG n 2 ). 

Definition. Let f n ,g n be two positive sequences. 

(1) We say that f n = @(g n ) if and only if there exist two positive reals a,f3 such that, for 
all sufficiently large n, ag n < f n < f3g n . 

fn 

(2) We say that f n ~ g n if and only i/hnin^oo — = 1. 

9n 
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5.1. Taylor series I. Let dfj 

zero vector) and S n = Yli<j Kj 

Proposition 8. OG n +i 2 



d°° 



+ h%j, h 



Proof. Let us recall that OG n 



Let e n 

J- 
16 

Sn 

16 



(h 



12 



(tiij)i<j (the sequence (h n ) converges to the 

n \2 joo 



'.31 



'12 



+ (h 



13 



"23^ "14- 



16 



With the help of the software "FGb" we put OG n+ {' 



in the form of a rational function in the unknowns (hfj) and we deduce the terms of degree at 
most two of its Taylor series: 



OG n 



+i 



16 



Sn 

1 + JL 
4 



+ 0{\\h n f) 



16 



1 

16 



+ 0(\\hn\\ z ) where: 



4 16 lV 12 



^ 4 ) 2 dg + (^ 3 



/?.. 



5.2. Taylor series II. Let us recall that dfl 
Now we prove the second key. 

1 2 \ 



12 



4 

"23 

• df| and df^ + df z 



2 df i )+0{\\h n \ 



□ 



> 4. 



Proposition 9. i) 5 n = 0(||/i n | 
ii) There exists k > smc/i i/iat e n > 

Proof. For i): From assertion ra) of proposition |2|, /i n satisfies an algebraic relation. We obtain 
with Maple the Taylor series of the preceding relation; we consider the terms of degree at most 

| 3 ) = where a' n = 0(\\h n \\ 2 ) 
0{\\K\\ 2 ). 

/134, then we transform the expression a' n 
^ — A)5 n + 0(\ \h n \ | 3 ). a n is a quadratic 
(/i™ 3 , hi^, h^z, /i24' ^34) ^ K 5 whose symmetric associated matrix E is defined by: 



-kS n - 



two: 4(4 - d°°)(4 - d^)(d^ + dg - 4)<5 n + < + 0(II0 3 ) = where of, = 0(11^,11' 



*12A-* - "13A"12 

negative quadratic form in the 

For ii): If we take h™ 2 
in a n that satisfies : a n ■ 
form in /it 



Thus S n 



/£n 



£ = 



£12 
£22 



£12 
2£i2 — £22 
£22 

* 



with < 



£11 
£12 

£14 = -dg a +4g^ 



"13 
-4(4 - d 

n un \ 
24' ' J 34J 

£l4 

£12 

£12 

£ll 

* 



n 14 



"23 



"24 



?§)(4-dS)(dS + d5 

G 



2£i 2 \ 
2£i2 

2£l2 

2£i2 
4£i 2 / 



32 - 12d^ + djjf - I2df 2 + d^ 



-Adf 2 + 16 



JOO JOO 

"12"13 



JOC 

"13 



-8rf 

4c? 



13 " 

X) 
13' 



'13 



By the same way we transform the expression e n in e' n such that e r , 
quadratic form in h' n whose matrix is: 



Odl^nH 3 ). 



E = 



with 





E\2 




E14 


2E 12 




E21 


2Er2 — E22 


E\2 


2E 12 
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E22 


E\2 


2E 12 


* 
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En 


2E\2 


V * 
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4£i2 


f E n 


— a 12 








1 E12 


— a 12 








| ^14 


— "12 


— ^00 

"13 






[ E22 


= 8- 


JOO 

"13 







\ 



un Ln _ in 
"24' "14 — "23i ' i 34 



-h 



satisfy 



Let -F be the plane of R defined by the relations: {/i™ 3 
i 7 is the nullspace of the matrices £ and -E 1 . 

Definition. t-Fe say t/iai /i n has an acceptable value if and only if the (d^ ) = (d, 
Proposition [| Hi). 

Assume h n has an acceptable value and a small norm; an easy computation proves that 
h' n G F if and only if 7^ is an isosceles tetrahedron. 

We know that if Tq is not isosceles then, for all n, the tetrahedron T n is not isosceles. Thus here 
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h' n ^ F and if h' n = u n + v n is the decomposition associated to M 5 = F © F- 1 -, where F 1 - is the 
orthogonal of F, then (u n ) tends to and, for all n, v n 0. Moreover h' n T 'Eh' n = v n T T,v n > 
and h' n Eh' n = v n T Ev n > 0. 

spectrum(E\ F± ) = {2dJ§, Sd 1 ^, 2dJ|} and the associated eigenvectors are [— 1, 0, 0, 1, 0] T , 
[1, 1, 1, 1, 2] r and [0, -1, 1, 0, 0} T . 

spectrum ( S , F± ) = { 2 (4 - df 2 ) (4 - dg ) , 8 (4 - dg ) (4 - dg ) , 2 (4 - df 2 ) (4 - dg ) } and the associated 
eigenvectors are also [—1,0, 0, 1, 0] T , [1, 1, 1, 1, 2] T , [0, —1, 1, 0, 0] T (E and £ are simultaneously 
diagonalizable) . 

JOO JOO JOO 

"13 "12 "14 



Let to = min 
M = max 



(4 - df§)(4 - dg) ' (4 - d°§)(4 - d?|) ' (4 - dg)(4 - df§) 

JOO JOO JOO 

"13 a 12 a 14 



(4 - df§)(4 - dg) (4 - df§)(4 - df2) (4 - dg)(4 - df§) 

e' (/&' ) 

Thus if /i n has an acceptable value then to < n n < M and 

e' n = e n + 0(\\h n \\ 3 ) > ma n = -4m(4 - dg)(4 - dg)(4 - d^)5 n + 0(||/i n | | 3 ). 

To finish the proof of ii), it remains to show that {e n and 8 n are @(||^n|| 2 )} (*)• 
Indeed (*) implies that for all sufficiently large n, e n > —k5 n , choosing a positive number 
k < 4p with p = min{dg(4 - dg), dg(4 - dg(4 - dg)}. □ 

3 

Remark, // i/ie /imii /orm is a regular tetrahedron then to = M = —. 

32 

ii ) It can be proved that p < — with equality if and only if the limit form is regular. 

In the following we prove the property (*). 
5.3. Improvement of certain estimates. 

Lemma 1. If (ijkl) G U = {(1234), (1324), (1423)} then dfydfo - (d°°) 2 = 0{\ \h n \ | 2 ). 

64 

Proof. Let L n = n = -=== = — . (L n ) tends to 1. 

("V ^12^34 + V ^13^24 + V ^14^23) 

We have (d°°) 2 = L„d$d&, d£d& - (d°°) 2 = (-?- - l)(d°°) 2 and, consequently, J- - 1 = 

VV 12 34 V 13 24 V 14 237 . Thus,ifwesetu n = v^^+ydV^+ydp^-8, 

we have to prove that u n = 0{\\h n \\ 2 ): 

Eameu y ^) 2 + ffi^- + fr&) + o(II M 2 ) - s 

+ 0(\\h n \\ 2 ) - 8 



u,. 




d9° 
fan _|_ /j n 

= = 0(| \h n \ | 2 ), by proposition ||. 

□ 

Lemma 2. //(ijifeQ G U then hfj + h^ = 0{\\h n \\ 2 ). 

Proof. d™d n kl - (d°°) 2 = dfiihtj + h kl ) + h t3 h kl = 0(\\h n \\ 2 ) by Lemma 1. □ 
Lemma 3. For all sufficiently large n there exists {ijkl) G IA such that 

\un _ un I > II 



Proof. By Lemma 2 there exists A > such that, for all n and (ijkl) E U, 

\hfj + h% t \ < A\\h n \\ 2 . Let e G ( '^[)- There exists 2V such that n > iV => \\h n \\ < e. Let 



n > N be a fixed integer. We may assume |/t" 2 | = supj <:? - Thus |/i™ 2 | > — 7=- > — 7=- and 



\ h h + h u\ < 6^4(^i 2 ) 2 and \h? 2 - h% A \ > 2\h% 2 \ - \hfa + h% A \ > 2\h^\ - QA(K{ 2 ) 2 . 

< \K[ 2 \ < \\h n \\ < e < then 6A(/i? 2 ) 2 < and \h% 2 - > |/# 2 I > TygllM- D 

Proposition 10. If (h n ) has acceptable values, then for all sufficiently large n, e n > A||fo n || 2 
inf iK j df- 

with A = — ; moreover —S n = Q(\\h n \\ ). 

6 

Proof. e n = (h% 2 - h^) 2 d^ 2 + (h% 3 - h%±) 2 df 3 + (7tf 4 - /i 23 ) 2 d^ and Lemma 3 give the first part. 
For the second part, if h n has an acceptable value then we know, by proposition ^, that: 

i) e' n = e n + 0(||/i n || 3 ). Thus, by the first part, if fj, < A then for all sufficiently large re, 

e'n > lA\ h n\\ 2 - 

ii) For all n, e' n (h' n ) < Ma n (h! n ). Then a n > J^\\hn\\ 2 (we can see a n as a function of h n ) for 
all sufficiently large n. 

Hi) -S n = — +0(\\h n \\ 3 ) with v = 4(4-dg)(4-d5§)(4-c^2). 

Thus if /^i < A, then for all sufficiently large n, — <5 n > \\h n \\ 2 . □ 

5.4. The main result in dimension three. 

Let r = max {MLl^i, ML^i , ML^i j £ ^ l} 

Theorem 3. XTie sequences (T2i)ieN an d (T^i+i^eN are we// defined and converge to two non 
planar isosceles tetrahedra that are symmetric with respect to O. Moreover the convergence is 
with at least geometric speed. 

Proof. OG n+1 2 = ^-^e n + 0(\\h n \\ s ) < ^ + ±k6 n + 0(| \h n \ | 3 ) ~ OG 2 (l - A 

(according to the preceding proposition). Finally, for all sufficiently large n, OG n+ \ < qOG n 

where q < ^jl — ^ = r. 

For all sufficiently large n, G n and G n+ i are close to O then close to the middle points of 
the segments L4 n ^4 n+ i] and L4 n+ i^4 n+2 ]; thus for all sufficiently large n, A n A n+2 < 3G n G n+ i. 
Let p G N*; A n A n+2p < 3 E^? -2 G fe G fe+1 < 3 E^?" 2 ^ + OG fc+1 ) < 6Efci?OG fc . 
The series E OG n converges, then (A 2n ) n is a Cauchy sequence; thus it converges to A 00 . By 
the same way (A 2n+ \) n converges to ^4'°°, the symmetric of A 00 with respect to O. Moreover 

A2nA°° < 6 JX 2 n OG fc < j^—OG 2 n for all sufficiently large n; therefore A 2n A°° = 0(q 2n ) 

and the sequences (T 2 i) and (7~2i+i) converge with at least geometric speed. □ 

5.5. About the limit form. We consider the system of barycentric coordinates, for M 3 , 
determined by the vertices Aq, Bq, Cq, -Do of Tq. Let(u,v,w,h) be four non all equal reals. 
{aAo + (3Bq + 7C0 + 5Dq ; ua + v(3 + + h8 = , a + /? + 7 + <5 = l}isa plane II. 

(u, v, w, h) are said to be the barycentric coordinates of II. 

We use the notations of 2.1. Tq is not isosceles; then {a'bc,b'ca,c'ab,a'b'c'} are valid coordi- 
nates and define a plane LT (Lemoine's plane) that does not intersect S. Let I±, I2 be the two 
points satisfying the following conditions: they are symmetric with respect to II and inverse 
with respect to S. They are said to be the isodynamic points of Tq. There exist only two inver- 
sions leaving invariant S and mapping Aq, Bq, Co, -Do on the vertices of an isosceles tetrahedron; 
ii,/2 are the centers of these inversions (see [[ ]] p. 184-186). The link with our paper is that 
the associated isosceles tetrahedra (that are not symmetric with respect to O) are isometric to 
the limits of the sequence {%) but these couples are distinct. 
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Remark. • The limit form is a regular tetrahedron if and only if the parameters of the non 

64 
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isosceles tetrahedron Tq satisfy ^12^34 = ^13^24 = ^14^23 < ^ that is Tq is an isodynamic 



tetrahedron. 

• According to numerical experiments, we conjecture that OG„+i ~ r x OG n 
(convergence with order one); in particular it can be observed the following facts: 

i) If the limit isosceles tetrahedron is almost flat, that is close to a rectangle, then the factor 
r is close to 1. 

ii) If the limit form is a regular tetrahedron then OG n +i ~ — x 0G n . 

6. Sequence of cyclic quadrilaterals 

6.1. Degenerate simplices. Now we work in the Euclidean plane and we have a look on the 
case where Tq is a cyclic quadrilateral. We see the following as a degenerate case of the preceding 
study. 

Theorem 4. We consider a convex cyclic quadrilateral Tq = {AqBqCqDq) with circumcircle C, 
the circle of center O and radius 1, and such that its vertices are pairwise distinct; we use the 
preceding iteration and produce a sequence of convex quadrilaterals (Ti)i. The sequences (T2i)i 
and (^2i+i)j are well defined and converge to rectangles that have same image, whose centroid 

is O and whose lengths of the edges are df^ = 4, = ^ 12 ^ — and dj| = 4 — d^. If the 

V dl3<i24 

- 21 

limit form is not a square then OG n +\ ~ — — OG n and the sequence converges with at least 

geometric speed; else the convergence is even faster. We study examples where the limit form is 
a square and the convergence is with order three. 

We keep the preceding notations. In particular the parameters {d%j)i<j refer to the square of 
the lengths of the edges (diagonals included) of T n . 

Proof. • The beginning of the proof is similar to that of Theorem |2[ We use Theorem [I] (note 
that inequality in) is useless because Pt(T n ) = 0); with ii) we prove that the sequences {dfjd^n 
are increasing and if i < j then the (dfj) n have a positive lower bound; thus a cluster point has 
pairwise distinct vertices. With i) we prove that a cluster point T is a rectangle. Assertion ii) 
of Theorem || is always valid and gives explicitly the (df-)i<j, the parameters of T. 
The quadrilaterals are convex, then a"^ = 4 and d^ A = 4 — d^ . Let fj, = • 

• Relations between the hu. 

Part I. There exist three algebraically independent relations between the 6 parameters (^)i<i- 
The computation of the terms of degree 1 of the Taylor series of the relation T(A n , B n , C n ,D n ) = 
gives h r { 2 + h% + h% 3 + - h% 4 - fey a = 0(\\h n \\ 2 ). 

Moreover the triangles (A n B n C n ), (B n ,C n , D n ) admit a circumscribed circle of radius 1; 
for example for the first triangle the relation is a n 2 b n 2 c n 2 = 16 x S n 2 (the value of the area 
S n of the triangle (A n B n C n ) is given in Proposition || i)). The computation of the terms of 
degree 1 of the Taylor series of these relations gives h™ 3 = 0(\\h n \\ 2 ) and h^ = 0(\ \h n \ \ 2 )), thus 
8 n = 0(\\h n \\ 2 ). 

Definition. We say that h n has an acceptable value if and only if the (dfj) = (d?? + hfj) satisfy 
the preceding three relations. 

Using the proof of Lemma ||, we show that h\ 2 + h 3A = 0(\\h n \\ 2 ) and h™ 4 + h% 3 = 0(\\h n \\ 2 ). 
Thus h1 2 2 + h1 2 = ®{\\h n \\ 2 ). 

Part II. The computation of the terms of degree at most two of the Taylor series of the relations 
associated to the triangles (A n B n C n ), (B n C n D n ) give: h™ 3 fi(4: — (jl) + (h^ + ^23) 2 = 0(||/i n || 3 ) 

— (h n — h n ) 2 

and /i>(4 - n) + (h% 3 + h% 4 ) 2 = 0{\\h n \ | 3 ) that is = ^4 _ ) + C(||^|| 3 ) and 

h n _ ~( h 12 + Ki) 2 , 0( n h n 3] 
h 2A~ m(4 _^ +G(\\h n \\). 



The relation T(A n , B n , C n , D n ) = gives by a similar way 
M/i - 4)(^ 2 + + h& + h& - - h%) + /i^ 4 2 + (4 - v)K n l2 2 = 0{\\h n \\ 3 ). 

j^n 2 hi 2 

Using the last three relations we can deduce: 5 n = ^- - + 0{\\h n \\ 3 ) and 

5 n = 0(\\h n \\ 2 ). Moreover e n = 4/^ 2 2 + 4(4 - n)K A 2 + 0(\\h n \\ 3 ). 
Finally OG n+1 2 = ( ^! + p£_) + 0(\\h n \\*) = { -^-OG n 2 + 0(\\h n \\ 3 ). 
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Cose 1: /U 7^ 2; then the limit form is not a square and OG n +i ~ — - — OG n . The convergence 
is with order one. 

3 3 

Case 2: fi = 2; then the limit form is a square and OG n+ i = 0(||/i n || 2) or OG n+ \ = 0(OG n 2). 
In fact the convergence is much faster than in the preceding case. For example if the polar angles 
of Aq, Bq, Co, -Do are {0, arccos(0. 923827833284), arccos(— 0.8), — arccos(0. 9)} then we obtain a 
quasi square after few iterations. 

Now we can conclude: in both cases the series ^ n OG n converges; thus we may reason as 
in Theorem |3] and we obtain that the sequences (T^ieN and (72i+i)ieN converge to rectangles 
that have same image with at least geometric speed. □ 

6.2. A particular case. The nullspace N of the symmetric matrix £ is the hyperplane defined 
by n^g = /124. Assume h n has an acceptable value and a small norm; then a geometric argument 
or an algebraic computation shows that h' n G N is equivalent to T n is an isosceles trapezoid. 

Proposition 11. Preserving the assumptions of theorem 4, we assume there exists k such that 
Tk is an isosceles trapezoid but not a rectangle. If moreover = 2 then OC n ,_|_i ~ OG n 3 . 

Proof. For all n > k, T n is an isosceles trapezoid that admits OGk as a line of symmetry. 
Therefore T n cannot turn around O; we know explicitly the limit form of T n , then T n converges 
to a rectangle that admits OGk as a line °f symmetry. 

We study the rate of convergence using an explicit calculation. We may assume that the 
line of symmetry is the axis of abscissas. The vertical edges A n D n and B n C n have a n ,b n as 

abscissas; then the abscissa of the centroid of T n is g n = - ^ — . An easy calculation gives: 

a n 2 b n + 2a n b n 2 + b n 3 - 4a n 

a n+l or,i or 2 , a ' 

a n z - 2a n b n - 3b n +4 



b n 2 a n + 2b n a n 2 + a n 3 - 4b n 



+ 1 - — r~2 



9n+l 



b n - 2b n a n - 3a n 2 + 4 
- b n tends to a: 

(a n + b n ){ai - 4a n 3 b n - 10a n 2 b n 2 + 16a n 2 - Aa n b n 3 + 166 n 2 - 16 + b n A ) 



d°° 

We know that a n + b n tends to and a n 2 tends to the known expression 



-32 



i) df^ 7^ 2. The limit form is not a square and we obtain g n +\ ~ ( 1 j^-j 9n (convergence 

with order one). If d^ > 2 then, for all sufficiently large n, O €]G n G n +i[. 

ii) = 2. The limit form is a square that is the parameters of Tq satisfy the relation (di2) 2 = 

^14^23- Let a n = —= + u n ,b n = — + v n where u n and v n tend to and (u n ,v n ) ^ (0,0). 

v 2 v2 

We calculate the Taylor series of the preceding relation and we consider the terms of degree 
at most 2; we obtain 16y/2(u n — v n ) = — 8(3u 2 + 3u 2 + 2u n v n ) + 0(||(n n ,f n )|| 3 ). Therefore 

u n -v n ~ ^=(&4 + 3-u 2 + 2u n v n ). u n -v n = 0(\\(u n ,v n )\\ 2 ) thus u n + v n = @(\\(u n ,v n )\\). 

w A A •! +ll + . + (a + fe)(-4(u n + t; ra ) 2 + 0(||(n ra ,^ ra )|| 3 ) -4(u n + v n ) 3 

We deduce easily the estimate g n +i ~ ~ 

32 32 

or <7 n +i ~ g n 3 (convergence with order three) what is astonishing. 

For example if ao = 0.955, 60 = 0.12237784429, then we obtain a quasi square after three 

iterations. □ 
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6.3. About the limit form. As in the case of a tetrahedron there exist two inversions leaving 
C invariant and mapping Aq, Bq, Cq, Dq on the vertices of a rectangle; the centers of these 
inversions are inverse with respect to C. The associated rectangles are isometric to the limit 
of the sequence (%). In particular the limit form is a square if and only if the parameters of 
the non rectangular quadrilateral Tq satisfy d^d^n = c?i4(i23 < 4 that is Tq is a harmonic (or 
isodynamic) quadrilateral. 

Remark, i) If the limit form is a square then the convergence seems to be with order three; 
more precisely we conjecture that OG n+ i ~ OG n 3 . 

ii) It can be observed a strange phenomenon: if d = 3 and the limit form is close to a flat 
tetrahedron then we have seen that r ~ 1. If d = 2 then r is in general far from 1; in particular 
if the limit form is close to a square then r is close to 0. 

An explanation is that, if d = 3, then G n tends to O in such a way that the angles between 

> 7T 

OG n and the edges of T n are close to —; thus if d = 2, then we study a quasi orthogonal 
projection of the preceding G n that converges faster than this G n . 



7. Solution of the case d = 2 
Now we consider a dynamical system of triangles. 

7.1. Some remarks. 

i) This problem has been solved in ||| by one of the authors. The proof in 0] is partially 
geometric; now we give a complete proof that is essentially algebraic. Moreover we give in 
Proposition 13 a much better estimate of OG n . 

ii) (j> is not one to one: indeed, if T\ is a generic triangle, then there exists two triangles Tq 
such that 4>(Tq) = T\. 

7.2. The parameters. Let a = BoCo,b = BqCq,c = CqAq; the parameters of Tq are: s = 
a 2 + b 2 + c 2 , t = a 2 b 2 + b 2 c 2 + c 2 a 2 , u = a 2 b 2 c 2 . 

s 2 

Let us recall that the circumradius of Tq is 1; then u = At — s 2 and t > — . Moreover u = 16 x S 2 
where S is the area of {AqBqCq). 

OGq 2 = 1 - - and < s < 9. Therefore s = 9^Gq = 0^ {AqBqCq) is an equilateral 
9 

triangle. 

s 2 — 3t = a 4 — a 2 (b 2 + c 2 ) + 6 4 + c 4 — b 2 c 2 ; this is a polynomial in a 2 with discriminant 

2 2 

-3(6 2 - c 2 ) 2 < 0. Thus t < ^- and t = ^- O (A BqCq) is equilateral. 

7.3. Inequalities. Here Aq,Bq,Cq can be on a line but are not all equal; thus s > 0,t > 0. 

s\,ti,ui are the parameters of T\; we obtain: 

s 2 (6t-s 2 ) s A t(9t-2s 2 ) „ „ , 

si = — ,h = ^ > u i = 4t i ~ s i wn ere 

D = —4s 3 + 18st — 108i + 27s 2 . D > because the numerator of s± is positive; moreover D is 
bounded. 

3sf9 — s)(4:t — s 2 ) 

• si — s = — > 0. Moreover si = s if and only if s = 9 or u = that is AqBqCq 

is equilateral or flat. 

9t(4t - s 2 )(9(— - t){s - 6) 2 + s 2 (9 - s)(s - 3)) 

• tx-t= ^ . if s > 3 then t x > t. 

2 

• ~H± = (^-^j . ui>u^s 3 >D^is = s 2 (9 - s) + 18 - t \ (s - 6) > 0. 

s 2 s 3 

If s > 6 then v > 0. Now we assume s < 6; t > — implies that v > — . 

- - 4 F ~ 2 

Therefore u\ > u. 
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7.4. Convergence of the triangles. We assume Tq is not an equilateral triangle. s n ,t n ,u n 
refer to the parameters of T n . 

Proposition 12. Let T be a cluster point of the bounded sequence (T n ) n ; then T is a non flat 
equilateral triangle. Moreover the lengths of the edges of T n converge to 

Proof. • The sequence (u n ) is increasing; thus it converges to u°° > 0. Let T be a cluster point 
of the sequence (T n ). The sequence (s n ) is increasing; thus it converges to s°° > 0. With a 
proof similar to that used in the theorem 3 we show that: 

(1) s 00 ,^ 00 are parameters of T; therefore T is not flat. 

(2) T is an equilateral triangle and s°° = 9, u°° = 27. 

• Moreover for all sufficiently large n, s n > 3 and the sequence (t n ) is increasing then convergent; 
finally (t n ) converge to 27, the corresponding parameter of T. The sequences (s n ), (t n ), (u n ) 
converge to 9, 27, 27; therefore (a n ), (b n ), (c n ), the lengths of the edges of T n , converge tO y/3. □ 

Let a n 2 = 3 + h n ,b n 2 = 3 + k n ,c n 2 = 3 + l n ,5 n = (h n ,k n ,l n ). \\.\\ refers to the euclidean 
norm. Now we show an important estimate. 

Lemma 4. h n + k n + l n ~ —{h n k n + k n l n + l n h n ) ~ — ll^n|| 2 when n tends to oo. 

3 6 

Proof. u = t- 4s 2 44> 3(h n + k n + l n ) = -(h n + k n + l n ) 2 + (h n k n + k n l n + l n h n ) 

-h n k n l n . Thus 3(h n + k n + l n ) = -(h n + k n + l n ) 2 + (h n k n + k n l n + l n h n ) + 0(||<5 n || 3 ) = 

+ 

Ill 

h n k n + k n l n + l n h n = -{h n + kn + l n ) 2 - -(h 2 + kl + l 2 n ) -(h 2 + k\ + l 2 n ). □ 

Proposition 13. If G n is the centroid of T n then OG n+ \ ~ OG n 2 when n tends to oo. Thus 
the sequence {OG n ) converges with order 2. 

Proof. OG r 2 = n — g n ■ OG n +\ 2 = 1 ^sing Maple and Lemma 1, we obtain the 

Taylor series of OG n+ \ 2 with the precision 0(||<5 n || 5 ): 

OG n+ i 2 = where 
— 81 z 

N n = 81(h n + k n + l n ) 2 + 18{h n + k n + l n )(2h n 2 + 2k n 2 + 2l n 2 + h n k n + kj n + UK) + 0(\ \5 n \ | 5 ) 

= 27 '(h n + k n + l n )(3(h n + k n + l n ) - 2 + lnhn)) + 0(\\5 n \\ 5 ) 

~ -81(h n + k n + In) 2 - Finally OG n+1 2 ~ ^ K + kn + l ^ = OG n \ □ 

ol 

Remark. • We can deduce that there exists A £ (0,1), that depends upon a,b,c, such that 
OG n ~ A 2 ". If To is close to a flat triangle then A is close to 1. Of course if Tq is close to an 
equilateral triangle then A is close to (see Tq as the result of a large number of iterations.) 

• The result obtained in OG n +i = 0(OG n 2 ) is weaker and does not give the preceding 
estimate of OG n - 

7.5. The main result in dimension two. 

Theorem 5. The sequences of triangles ('?2i)iGN and ( / 72i+i)ieN ore well defined and converge 
with at least quadratic speed to two equilateral triangles that are symmetric with respect to O. 

Proof. For all sufficiently large n, G n and G n +i are close to O then close to the middle points of 
the segments L4 n A n+1 ] and [ J 4 n+ i J A n+ 2]; thus for all sufficiently large n, A n A n+2 < 3G n G n +i- 
Let peN*; AnA n+2p < 3 Ek= n P ~ 2 G k G k+1 < 3 El= n P " 2 (OG k + OG k+1 ) < 6 OG k . The 

series Yl OG n converges, then (A2 n )n is a Cauchy sequence; thus it converges to A 00 . By the 
same way (^2n+i)n converges to ^4'°°, the symmetric of A°° with respect to O. 

Moreover A 2n A cc > < 6 EfcL 2 n 0G k < 12 x OG 2n for all sufficiently large n; therefore A 2n A co = 
0(X 2 ") and the sequences (T 2 i) and (T 2 i + ±) converge with at least quadratic speed. □ 
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8. Conclusion 



We mention some questions we can ask about these dynamical systems. 

i) Can one generalize our results concerning the (i-simplices for d > 3? We note that the 
complexity of the computations increases very quickly 

ii) What occurs if Tq is a degenerate <i-simplex that admits a circumsphere in IR e with e < d? 
For example, Tq could be a convex polyhedron that admits a circumsphere in M 3 . 

iii) More generally what occurs if we replace the centroid of % with some barycenter of the 
vertices of %! 
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